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Abstract. We review the 3+1 - split which serves to put Einstein's equations into the 
1 form of a dynamical system with constraints. We then discuss the constraint equations 

under the simplifying assumption of time-symmetry. Multi-Black-Hole data are pre- 
, ^ , sented and more explicitly described in the case of two holes. The effect of different 

Q 1 topologies is emphasized. 

<' 

. Notation. Space-time is a manifold M with Lorentzian metric g of signature 

CN | (—,+,+,+). Greek indices are G {0,1,2,3} and latin indices are e {1,2,3}. 

Indices from the beginning of the alphabet, like a,/3, . . . and a, b, . . ., refer to 
orthonormal frames and indices from the middle, like A, fj,, . . . and I, m, . . . to co- 
ordinate frames. The symbol o denotes the composition of maps. The relation := 
1/-) | (=:) defines the left (right) hand side. The torsion and curvature tensors for the 

connection V are defined by T(X, Y) := V X Y - Vyl- [X, Y] &ndR{X,Y)Z := 
VxVyZ — VyV^Z — V[x,Y]Z respectively. The covariant components of the 
Riemann and Ricci tensors are defined by R a /3-yS '■= ff(e a , i?(e 7 , es)ep) and 
Ra-y ■— g /3A Rai3j8 respectively. 
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In this article we discuss the vacuum Einstein equations 
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7" 1 The 3+1 - Split 

> 

x. 

which form a system of ten quasi-linear second order differential equations for 
the ten functions g^ v . However, the four equations G M o = do not involve the 
second time derivatives and hence constrain the set of initial data. To see this, 
recall that the twice contracted second Bianchi identity gives V^G^" = 0, or 
expanded 

dvG au = -d k G kv - r^ x G Xv - r^ x G" x . (2) 

Since the right hand side contains at most second time derivatives the assertion 
follows. The ten Einstein equations therefore split into four constraints and six 
evolution equations G lk = 0. That four equations constrain the initial data rather 
than guiding the evolution results in four dynamically undetermined functions 

* To appear in Black Holes: Theory and Observation, edited by F. Hehl, C. Kiefer and 
R. Metzler, Springer Verlag 1998. 
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among the ten g^ u . The task is to parameterize the g^ v in such a way that four 
dynamically undetermined functions can be cleanly separated from the other six. 
How this can be done via the 3+1 split is explained below. The four dynamically 
undetermined quantities will be the famous lapse (one function) and shift (three 
functions). It follows directly from (2) that the constraints will be preserved 
under this evolution. 

The splitting of the Einstein equations will be formulated in a geometric fash- 
ion. We initially think of (M, g) as given and satisfying the Einstein equations. 
Then we write down the evolution law for the intrinsic and extrinsic geometry of 
a spacelike 3- manifold £ as it moves through M. Together with the constraints 
they are equivalent to all Einstein equations. Finally this procedure is turned up- 
side down by taking the evolution equations for £'s geometry as starting point. 
Only after their integration can we construct the ambient space-time. 

1.1 3+1 Split Geometry 

The topology of space-time (or the portion thereof) which we want to decompose 
into space and time must be a product M = £ x R. We foliate M by a one- 
parameter family of embeddings e t : £ — > M , t ? R. For fixed t the image of 
et in M is called £ t , or the i'th leaf of the foliation. All leaves are assumed 
spacelike. Hence there is a normalized timelike vector field, n, normal to all 
leaves. We choose one of the two possible orientations and thereby introduce the 
notions of future and past: A timelike vector X is future pointing iff g(X, n) < 
(recall signature convention). The tangent-bundle T{M) can now be split into 
the orthogonal sum of the subbundle of spacelike vectors, S(M), and the normal 
bundle, N(M). The associated projection maps are given by 



which can be naturally continued to the cotangent bundle by setting S*(u) := 
oj o S and then factorwise on tensor products and linearly on the whole tensor- 
bundle. Thus we obtain a split of the whole tensor bundle, where from now on 
the projection maps are simply called S and l\l for all tensors. Tensors in the 
image of S are called spatial. It is easy to verify that 



where n b := g{n 1 •). Note that the restriction h t of h to T(£ t ) is just the induced 
Ricmannian metric on £ t . Identifying for the moment £ and £ t via this 
leads to h t — e* t g (Exercise). For what follows it is however crucial to regard 
spatial tensors as tensors over M and not over £. Otherwise covariant (or Lie-) 
derivatives in directions off £ would not make sense. 
If X, Y are any spatial vector fields we can write 



S : T{M) — > S(M), X^X + ng(n,X) 
N : T(M) — > N(M), X^-ng(n,X), 



(3) 
(4) 



h := Sg = g + n b <g> n b , 



(5) 



V X Y = SV X Y + NV X Y = D X Y + n K(X, Y), 



(6) 
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where we defined the spatial covariant derivative, D, and the extrinsic curvature, 
K, by 

D x :=SoV x , (7) 
K(X, Y) := -g(V x Y, n) = -g(V Y X, n) = g(V x n, Y). (8) 

The second equality in (8) - and hence the symmetry of K - follows from the 
vanishing torsion of V and the fact that [X, Y] is spatial. It is easy to prove that 
K is indeed a tensor and that D defines a connection on the tangent bundle of 
each leaf S t . Extension via the Leibnitz rule leads to a unique connection on the 
bundle of spatial tensors, which can be directly defined by (7) with the extended 
meaning of S described above. In fact it is just the Levi-Civita connection com- 
patible with the metric h. To see this, we compute Dxh = SVx{g + n b <8m b ) = 0, 
since S7 x g = = Sn b , so that D is compatible with h. Vanishing torsion is 
also immediate: D X Y - D Y X - [X,Y] = S(V X Y - V Y X - [X,Y]) = 0, by 
[X, Y] = S[X, Y] and the vanishing torsion of V. 

Let {e , ei, e2, 63} be an orthonormal frame adapted to the foliation, i.e. 
e = n, and {e°, e 1 , e 2 , e 3 } its dual. Then from (5) with n b — e° we have 

3 

g = -e ®e° + h= -e° ® e° + ^ e a ® e Q . (9) 

a=l 

The family of embeddings t ^ e t defines a vector field, d/dt —: dt, which is 
easily characterized by its action on any smooth function /: 

dtf:=^. foe t . (10) 
at t=o 

This vector field can be decomposed into normal and tangential components 

d t = an + f3 = ae +f3 a e a , (11) 

with uniquely defined function a and spatial vector field (3. They are called the 
lapse (function) and shift (vector field) respectively. 

Let now be an adapted local coordinate system on M so that x° = t 

and hence spatial fields dk '■= d/dx k . The flow lines of dt are then the lines 
of constant spatial coordinates x k . Hence (a, (3) are interpreted as normal and 
tangential components of the 4-velocity - measured in units of t — with which 
the points of constant spatial coordinates move. To express the metric g in terms 
of these coordinates we use an obvious matrix notation and write 

^(;: ; )(::) : 

(e° e«) = (dt dx k )(l (13) 

Introducing (13) into (9) yields the 3+1 split form of the metric g: 

g = -a 2 dt®dt + h ik (dx l + (3 l dt) ® (dx k + f3 k dt), (14) 
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where hik = h(di, dk) = J2 a AfA^ and /3 % Af = (3 a . For the measure 4-form one 
easily obtains e° A e 1 A e 2 A e 3 = ctyJdct{h ik }d A x. 

In the ambient space-time the notion of time-derivative of spatial tensors is 
introduced via the Lie-derivative along the time flow generated by d t - But in 
order to render this an operation within the space of spatial tensor fields we 
must include a spatial projection. Using (11) we define the "doting" by 

h := SLg t h = aL n h + SLph, (15) 

where we also used L an h = aL n h and that L n h is already spatial. This is true for 
any covariant spatial tensor and any smooth function a. To prove this, we first 
remark that by Leibnitz' rule it suffices to prove it for a general spatial 1-form lo. 
The first assertion now follows from L an uj = (i an od+doi an )u) — ai n du) = aL n cj. 
The second statement follows from the general formula i n °L v — i[ n _ v ] +L v oi n , 
showing that for v = n the left hand side annihilates any spatial tensor field. 
This identity also shows why we need the projector in the second expression 
on the right hand side of (15), since for v = [3 it shows that we would need 
[n, ft] oc n for Lph to be spatial. But this is generally false, as one easily shows 
that [n,(3 k d k ] ocnO d t ((3 k ) = 0. 

We proceed by showing that L n h is just twice the extrinsic curvature: 

K = \L n h. (16) 

To prove this relation, we take any spatial vector fields X, Y and compute: 
L n h(X, Y) = V„(/i(X, Y))-h([n, X],Y)-h(X, [n, Y]) = h{V x n, Y)+h{X, V y n) = 
2K(X,Y), where we used the metricity of V, {V„h){X,Y) = (V„g)(X,Y) = 0, 
and its vanishing torsion. Hence we arrive at 

K = (h - SL (j h) . (17) 

The projected Lie-derivative can be expressed in terms of the spatial covariant 
derivative in the usual way. In components with respect to a spatial coordinate 
frame this reads (SLph)ik — Diflk + Dkfti- 



1.2 Constraints and Equations of Motion 

Using the splitting formula (6) for the connection V in terms of D and K we 
can derive the so-called Gauss-Codazzi and Codazzi-Mainardi equations by a 
straightforward manipulation. In components with respect to {e a } and with 
i?( 3 ) denoting the curvature of D, they read respectively: 

Rabcd = R^abcd + K ac Kbd — K a dKb c , (18) 

Roabc - D c K ab - D b K ac . (19) 

From here it is easy to write down the constraints by noting that in orthonormal 
frames one has J^abRabab = R + 2i? o = 2Gqo, i-c. the 00 component of the 
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Einstein tensor just depends on the spatial components of V's curvature. In 
fact, it is the sum of the spatial sectional curvatures of V. Further, Got = 
Rob = ^2 a Roaba- Hence we have the constraints, now written in components 
with respect a coordinate frame, 

G^nV = i(i?( 3 ) - K ik K ik + (K]f) = 0, (20) 
G^ = D k (K ik -h ik Ki)=0. (21) 

To obtain the dynamical equations one starts again from the defining equa- 
tion of the curvature and manipulates the expression for i?oaOb- Observing that 
V n (g{e a ,V eb n)) = L n K ab one arrives at 

Roaob = -L n K ab + K ac K b + a a a b + D a a b , (22) 

where a := V„n. Note also that a b = L n n b (Exercise: Prove it). Despite appear- 
ance, the last term in (22) is also symmetric. 1 Now, R ab — —RoaOb + Reach, 
so that with (18) we have 

Rab = R {3) ab + L n K ab + K ab K c c - 2K ac K c b - a a a b - D a a b . (23) 

This is almost the evolution equation we wish to obtain. As in (15) we have 
K = aL n K + SLpK, and since we want to write down the final equation in 
a coordinate basis, we can simplify the terms involving a by noting that at = 
{L n n h ){di) = n b ([d,, ±(dt - 13)}) = did/a. Hence 2 

K lk = a(2KijKl - K ik K] + R lk - R^ lk ) + L K ik + D t D k a, (24) 

where in the vacuum case we consider here one sets Ri k = 0. Note that in a coor- 
dinate frame dotting just means taking the partial derivative of the components, 
i.e., L dt hij = dhij/dt. 

The dynamical formulation is now complete. The constraints are given by eqs. 
(20) (21) and the six evolution equations of second order are written as twelve 
equations of first order, given by (15) and (24). The six dynamical components of 
g are the hij , whereas there are no evolution equations for the four functions a, [3. 
The initial value problem thus takes the following form: 1.) choose a 3- manifold 
E with local coordinates {x*}, 2.) find a Riemannian metric hij and a symmetric 
covariant tensor field Kij on £ which satisfy (20)(21), 3.) choose any convenient 

1 This is due to n being hypersurface-orthogonal. To see this, we first note the identity 
(L„ — T;i„od)dn An = da An (Exercise: Prove it). Now, hypersurface-orthogonality 
of n dn b A n b = da A n b = <4> Sda b = <S> D [a a b] = 0. 

2 Be aware that some authors define the extrinsic curvature with opposite sign, for 
example E. Seidel in his lecture. Hence the discrepancy between our eqns. (17) (24) 
with his (5) (6) respectively. In our convention, which agrees with Hawking & Ellis, a 
positive K\ implies volume expansion under deformations in normal direction. Also 
note that "dotting" does not commute with index raising. Hence notations like K l ° 
are ambiguous. For example, denoting the index raising operation by a superscript 
t), (16) immediately gives (L„(/r 8 ) - (L n K) t y k = -±K)K ]k . 
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functions a(t, x k ), (3 l (t, x k ), 4.) evolve hij and Kij via (15)(24) by using the 
choices made in the previous step, 5.) take the solution curve hij(t,x k ) and the 
functions from 3.) to construct the space-time metric according to (14), where 
a; = t. The g so constructed solves Einstein's equations. An important theorem 
guarantees that for suitably specified data a maximal evolution (M,g) exists 
which is unique up to diffeomorphisms (Choquet-Bruhat and Geroch 1969). See 
also Choquet-Bruhat and York (1980) for a review and further references. 

Regarding step 2.), we remark that all topologies E allow some initial data, 
i.e., there are no topological obstructions to (20)(21) (Witt 1986). This might 
change if geometrically special data are sought (see below). To illustrate step 
3.), we mention the so-called maximal slicing condition on a. To derive it, we 
compute L n (h ij Kij) = 2/\ '-'A , ; • /''••'/.„ /\, , = -R 00 - K i:j + Aa/a, where 
A = D l Di and where we used i?oo = So -^ooOa and (22) to replace L n K. Hence 

L dt (Kj) = (A- K^Kij - Ro a )a + L (Kj), (25) 

where we left in i? o for generality. Note that the strong energy condition implies 
-Roo > through the Einstein equations. The hypersurface E C M is called 
maximal if the trace of its extrinsic curvature - the so-called mean curvature - 
is zero, i.e., K\ — 0. This is equivalent to E being a stationary point of all 
the 3-dimensional volume functionals for domains in E and fixed boundaries. 3 
(Exercise: Prove this using (16).) Now, given a maximal slice E C M, (25) gives 
the following simple condition on a if the evolution is to preserve maximality: 
Oa = with elliptic operator O = A — K v — i?oo- (Exercise: Assuming the 
strong energy condition, prove that any smooth function a in the kernel of O 
cannot have a positive local maximum or negative local minimum on E.) In the 
vacuum case one can use (20) to write O = A — R( 3 \ i.e., purely in terms of the 
intrinsic geometry of S, where clearly > 0. 

The maximal slicing condition plays an important role in numerical evolution 
schemes, since - by definition - the evolving maximal slices S t C M approach 
slowest the regions of strongest spatial compression. In this sense they have 
the tendency to avoid singularities. For further information see section 2.3 of 
E. Seidel's lecture. Finally we note that since not all topologies E allow for 
metrics with R^ > 0, there exist topological obstructions to maximal initial 
data sets (Witt 1986). 

2 Time-Symmetric Initial Data 

Suppose a hypersurface E C M has vanishing extrinsic curvature, K = 0. From 
(6) we then have VxY = DxY for all vector fields X, Y tangent to E. In 
particular, if 7 : / — > E is a curve with tangent vector field 7' over 7, then 

3 The standard terminology is that such stationary points are called "maximal" if the 
ambient geometry is Lorentzian and "minimal" if it is Riemannian, irrespectively of 
whether they really are true maxima or minima respectively. True extrema are called 
stable maximal (minimal) surfaces. 
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V 7 <7' = -D 7 '7' and 7 is a geodesic in S iff it is a geodesic in M. Submanifolds 
for which this is true are called totally geodesic. This is a stronger condition 
than maximality. In general, constant mean curvature data play an important 
role in the solution theory for the constraints (see York 1973, Murchadha and 
York 1974). Here we shall vastly shortcut the general procedure by imposing 
the condition that E be totally geodesic. One can then show that the maximal 
development, M, from these data allows an isometry fixing E pointwise and 
exchanging the two components of M — E. Hence such data are called time 
symmetric. For such cases the constraints reduce to the simple condition that 
(E,h) has vanishing Ricci-scalar: 

R^(h) = 0, (26) 

where for later convenience we explicitly indicated the metric as argument of 
i?( 3 ) . A general idea for solving (26) is to prescribe h up to an overall conformal 
factor <P, and let (26) determine the latter. So setting h = <P 4 h', with fourth 
power just for convenience, we have by the conformal transformation law for the 
Ricci-scalar 

R(V($ 4 h') = -8$- 5 (A h , - ii?( 3 ) (//))<£ =: -8$- 5 C h ,$ = 0, (27) 

where ZW is the Laplacian for the metric h'. We are interested in C 2 solutions 
satisfying <P > and where (E, h) has no boundaries at finite distances, i.e. 
E should be topologically complete in the metric topology defined by the dis- 
tance function induced by h. The last condition is equivalent to (E, h) being 
geodesically complete (theorem of Hopf-Rinow-DeRahm, see e.g. Spivak 1979). 
In addition, we shall only be interested in manifolds whose ends are asymp- 
totically flat. Allowing the manifold E to have more ends or to be otherwise 
topologically more complicated allows for a greater variety of solutions. Note 
that to each of n asymptotically flat ends there corresponds an ADM-mass of 
which n — 1 are independent (see below). 

Brill Waves. One may ask whether simple asymptotically flat solutions to 
Ch'^ = exist on E = R 3 . There are no (regular!) black-hole solutions with 
this simple topology, but there are solutions representing localized gravitational 
waves of non-zero total ADM energy (Araki 1959). In the axisymmetric case they 
were investigated in detail by Brill (1959). Solutions of this kind are collectively 
called "Brill waves" . One takes (from now on in the usual shorthand suppressing 
the ®) 

ti = aq>(\q(z,p))(dz 2 + dp 2 ) + p 2 dip 2 , (28) 

where the profile-function q must for r — > 00 fall off like r~ 2 and like r~ 3 in its 
first derivatives in order for h to turn out asymptotically flat, q characterizes the 
geometry in the meridial cross section (zp-planc) of the toroidal gravitational 
wave. Regularity on the axis also requires q and d p q to vanish for p = 0. The 
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parameter A 6 R + is sometimes introduced to independently parameterize the 
overall amplitude. Equation (26) for <&(z,p) takes the particularly simple form 

(At + \\AMq) <P = 0, (29) 

where Af is the flat Laplacian and A^ = d 2 /dz 2 + d 2 /dp 2 . Given q, everywhere 
positive solutions for <P exist provided A is below some critical value depending 
on the choice q (Araki 1959). To see uniqueness, assume the existence of two 
solutions <?i and $ 2 and set hi = @fh', i = 1,2. Then <P 3 := <Pi/<£ 2 is also 
C 2 , positive and tends to 1 at infinity. But (27) immediately implies Ch 2 &3 = 
-g#!-R (3) (/ii) = 0, and since also i? (3) (/i 2 ) = this is equivalent to A h2 <P 3 = 0. 
Hence $3 = 1 due to the fact that the only bounded harmonic functions are the 
constant ones. 

3 Black-Hole Data 

A substantial variety of time-symmetric black-hole data can already be obtained 
by solving (27) when h' is flat, i.e., where the 3-metric, h, on the spatial slice 
at the moment of time-symmetry is conformally fiat. One can obtain manifolds 
with any number of asymptotically flat ends, and then reduce this number by a 
process which is best described by calling it "plumbing" (see below). We shall 
devote the rest of this paper to the description of such solutions and techniques. 
Note that for flat h! we are left with the simple harmonic equation involving 
only the flat Laplacian: 

A { <P = 0. (30) 

In general it is difficult to infer from given initial data whether they cor- 
respond to a spacetime with black holes, i.e. with event horizons. However, in 
the examples to follow it is easy to see that that there will be apparent hori- 
zons, since for time symmetric data apparent horizons correspond precisely to 
minimal surfaces S C E 4 . Proposition 9.2.8 of Hawking and Ellis (1973) now 
implies the existence of an event horizon whose intersection with E is on, or 
outside, the outermost apparent horizon for any regular predictable spacetime 
that develops from data satisfying the strong energy condition. Concerning the 
topology of apparent horizons we remark the following: Using the formula for 
the second variation of the area functional and the theorem of Gaufi -Bonnet, 
one shows that for ambient metrics with non-negative Ricci scalar any connected 
component of an orientable stable minimal surface of finite volume must be a 

4 The condition on S being an apparent horizon is that the congruences of outgoing 
null rays from S must have zero divergence. Analytically this translates into tr2(/t) = 
i(tr(K) — K(y, v) where k, v are respectively the extrinsic curvature and normal of 
S in S. The upper sign is valid for past apparent horizons, and the lower one for 
future apparent horizons. tr2 is the 2-dimensional trace using the induced metric of 
S and tr the 3-dimensional trace using h. For time-symmetric initial data (K — 0) 
this condition states that k is traceless and hence S minimal in E. 
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topological 2-sphere (Gibbons 1972). Allowing also for non-orientable apparent 
horizons, one deduces from this that for metrics h of £ with R^(h) > a 
connected component of an apparent horizon is either S 2 or RP 2 , the latter 
being the (non-orientable) 2-dimensional real projective space. If £ is orientable 
RP 2 C £ is one-sided, as in the example below. 

3.1 Schwarzschild Data 

We start by noting that the most general non-trivial solution of (30) on £ = 
R — {0} is given by ^(x) — 1 + ^ with r — ||x|| and to e R + . We cannot have 
any higher multipole moments because then necessarily has zeros on £. Just 
removing # _1 (0) from £ docs not work since these points are at finite distance 
so that the resulting space would not be (geodesically) complete. This is also 
the reason why to must be positive. Hence we obtain for the metric h in polar 
coordinates 

h= (l + g) 4 (dr 2 +rW), (31) 

with dQ 2 — d9 2 + sin 2 6 dtp 2 . Now, it is easy to verify that the following two 
diffcomorphisms, / and /, of £ are involutive (i.e., square to the identity) isome- 
tries: 

i(r,e,<p)~(^:A<py (32) 

I(r,6,<p) ■= (^T^-^ + t) ■ (33) 

The map / is called an inversion on the sphere r = to/2, whereas / is that 
inversion plus an additional antipodal map on the spheres of constant r. We shall 
sometimes refer to them as inversions of the first and second kind respectively. / 
has no fixed points while / fixes each point of the sphere S — {x | r = y } (which, 
as set, is also left invariant by /). As fixed point set of an isometry S must be 
totally geodesic 5 , hence minimal and therefore an apparent horizon. Its surface 
area is A = 167tto 2 , and it separates the two isometric regions r > to/2 and 
r < to/2. The metric (31) corresponds to the spatial part of the Schwarzschild 
metric of mass to in isotropic coordinates, which cover both asymptotically flat 
regions (I and III) on the Kruskal manifold. Using this isotropic form, one can 
read off a = (1 — m/2r)/(l + to / It ) , (3 = and verify that with this choice the 
static form of (24) with K = is satisfied (Exercise). 

The manifold £ has two isometric ends and we can get rid of one by suitably 
identifications. For this we take the quotient £ of £ with respect to the free 
action of /. The freeness guarantees that the quotient will be a manifold, and, 
by being an isometry, the metric descends to a smooth metric on the quotient. £ 

5 Proof: Consider the unique geodesic 7 starting on and tangentially to S. It cannot 
leave S since if it would, its image under / would be a different geodesic with the 
same initial conditions, which contradicts the uniqueness theorem for ODE's. 
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can be pictured by cutting E along S, throwing away one piece, and identifying 
opposite points on the inner boundary S on the retained piece. Hence topolog- 
ically E is the real projective space, RP 3 , minus a point. The projection of S 
into E is a totally geodesic, one-sided (i.e. non-orientable) surface diffeomorphic 
to RP 2 . E is orientable, smooth, complete and with one end which is isometric 
to, and hence has the same ADM mass as, cither end in E. This demonstrates 
how the introduction of more ends or other topological features makes it possible 
to define non-trivial black-hole data. One may also combine Brill waves with a 
black hole to model a single distorted black hole. This is further discussed in 
section 2.2 of E. Seidel's lecture. 

Multi-Schwarzschild Data. Taking E = R 3 — {ci, ■ ■ ■ ,c n } the generalization 
of (31) is easily obtained with n poles of strengths a* G R + at "positions" 

n 

*(x) = l + £* (34) 
»=i n 

where := ||x — c,||. For each i we can introduce inverted polar coordinates 
r 'i = a i/ r i to probe the region — > by letting r\ — > oo. Doing this shows 
that the metric is asymptotically of the form (31) with certain mass parameters 
m = rrii given below. The same is true for the region r — > oo with mass M. 
Hence one obtains n + 1 asymptotically flat ends. The internal masses and the 
overall mass are given by (r^ := ||cj — Cj||) 

rrii = 2oi(l + Xi), where Xi '■= ^2 ~> and M = 2^^. (35) 
In terms of the parameters a^r^ the binding energy takes the simple form 

n n n 

AM := M - ^ rm = -2 ^ a lXl - -2 £ ^ ^ < 0. (36) 

Note that there are as many independent masses as there are generators of the 
second homology group of S. These generators may be represented by stable 
minimal surfaces associated to each internal end. Their surfaces areas clearly 
satisfy A, > 167r(2ai) 2 , since the right hand side represents the minimal area 
in the strictly smaller metric (31) for just one hole with parameter m = 2a». 
But there is also an upper bound for the area, given by the recently proven 
Riemannian Penrose inequality 6 (Huisken and Ilmancn 1997), which in our con- 
text reads A < 167r(2ai) 2 (l + Xi) 2 ■ Assuming the existence of an event horizon 

6 The proof of Huisken and Ilmanen (1997) applies to all asymptotically flat Rieman- 
nian 3-Manifolds whose Ricci scalar satisfies R > 0. They prove that the area A of 
the outermost stable minimal surface bounding an end and the ADM mass m of that 
end satisfy A < l&Tim 2 . It implies the positive mass theorem for data with R > 0. 
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(see above), the Area Theorem (see my other contribution to this volume) im- 
plies that the area of the hole in the i'th end cannot evolve below Ai, which, 
using the first inequality above, implies in particular that the energy which is 
bound in the final hole is greater than 2a^ = + Xi)- The difference of the 

(conserved) ADM mass to the mass of the final hole is therefore bounded 
above by mjXi/(l + Xi)- I n other words, the fraction of energy being radiated is 
less than Xi/(1 + Xi)- This still allows for total conversion into radiation if one 
chooses Xi 00 • 

It would be of course more interesting to express (36) in terms of physical 
variables, like the individual masses rrii, and more geometrically defined distance 
functions than ry, like e.g. the proper geodesic distance of the minimal surfaces 
in the i'th and j'th throat. Note that for small mass-to-separation ratios we 
may in a first approximation replace cij by \rrii and by the geodesic distance 
of the i'th and j'th apparent horizons and get the familiar Newtonian formula. 
But there will be corrections the precise form of which depend on ones definition 
of "distance between two holes". Whereas here mass is unambiguously defined 
for each hole (by ADM), there is no natural definition of distance. Perhaps 
the easiest intrinsically defined distance is the one given above. For the multi- 
Schwarzschild manifold it has the disadvantage that the minimal surfaces are 
not easy to locate analytically and one has to resort to numerical methods (see 
Brill and Lindquist 1963 for early attempts). 

The location of minimal surfaces is interesting for a variety of reasons. It 
somewhat simplifies in the case of just two holes, which is automatically ax- 
isymmetric. Then the variational principle for the minimal surfaces reduces to a 
geodesic principle for curves in the zp- half-plane (cylindrical coordinates). The 
appropriately parameterized solution curves just describe a motion of a point 
particle in the potential — \p 2 $ & (Cadcz 1974). However, general analytic so- 
lution still do not exist. Numerical studies by Bishop (1982) for equal masses 
(ai = <i2 =: a) show the very interesting behaviour above the critical value 
a/r\2 — 1/1.53, where two more minimal surfaces appear, each of which en- 
closing the previous two. Initially they coincide, but for increasing a/r\2 they 
separate with the inner one rapidly increasing in area whereas the outermost 
staying almost constant. See also Gibbons (1984) for a related discussion. 

For the data discussed below the difficulty of determining location and size of 
minimal surfaces is absent, but somewhat as trade-off the concept of individual 
mass now becomes slightly more problematic. 

Different Topologies for Multi-Hole Data. There are other generalizations 
of the single hole case. The ones we discuss now will preserve the existence 
of involutive isometries like (32-33), but now for each apparent horizon. The 
manifolds they exist on have two or even just one end. The construction is 
somewhat involved (Lindquist 1963) and uses the method of images to construct 
solutions to (30). This method was introduced by Misner (1963) for the time 
symmetric case and later generalized to more general situations (e.g. Bowen and 
York 1980, Bowen 1984). (There is also a recent alternative proposal by Brandt 
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and Brugmann (1997).) For the general understanding it will be sufficient to 
explain the construction for just two holes. Note that the ADM definition of mass 
cannot be applied to the individual hole if it does not have an asymptotically 
flat end associated to it. But there exist alternative proposals for mass due 
to Lindquist (1963) and Penrose (1982) which can be employed here. (See also 
the general review by Penrose (1984).) But it should be pointed out that these 
definitions do not always apply in more general situations. For example, for the 
applicability of Penrose's mass definition within time-symmetric hypersurfaces 
the metric of this hypersurface must be conformally flat (Tod 1983, Beig 1991). 



3.2 Two Hole Data 

Just as in electrostatics, we shall use the method of images to construct special 
solutions to (30). This is done by placing image masses in an auxiliary, fictitious 
space so as to enforce special properties of <P. The properties which will be en- 
forced here are such that the inversions (32) (33) on 2-spheres become isometries. 

We start by drawing two 2-sphercs Si := S(ai, c,), i = 1,2, with radii a,i and 
centered at Cj. The spheres are non-intersecting and outside each other, so that 
ri2 > a\ + a 2 . On R 3 — {c^} we have the diffeomorphisms li and /j, which in 
polar coordinates at take the forms (32) and (33) respectively. These induce 
involutions on the space of functions, defined by 

Ji(f)--=-f°Ii and Mfl-^foii (37) 

n n 

respectively, where / is any function. The crucial property of these maps is 

Af o Ji = (fli/r^) 4 Ji o A{ and A f o J { = (ai/n) 4 J { o A f , (38) 

which in particular implies that the image of a harmonic function will again 
be harmonic, although with different singularity structure. The image of the 
constant function, / = 1, under either of these maps is just /' = ai/ri, i.e., the 
pole of strength a, at Cj. Moreover, given the unit pole /(x) = l/||x — d|| at d 
outside Si, then its image under Jj is 

Ji(/) = i^^diiiix-^d)!!' (39) 

and correspondingly for Jj. It represents a pole of strength a,/||cj — d|| < 1 at 
the image point Ii(d) (resp. Ii(d)). 

Writing down the metric h = <P 4 dsj in polar coordinates centered at a 
one easily verifies that 7j (7j) is an isometry of h if <P is invariant under Ji 
(Ji). The construction of such an invariant <P is by brute force: One averages 
the function 4>o = 1 over the free product of the groups generated by J\ , Ji 
{Ji, J 2)- The elements of this free-product-group are strings of alternating Ji's 
and J2S, where for each string length n > 1 there are the two different elements 
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Ji o J 2 o Ji ■ ■ ■ and J 2 o J 1 o J 2 - ■ ■. By definition, the string of length is the 
identity element. Hence one sets 

N 

*iv:=l + ££jii -"°Ji„(*o), (40) 

n=l 

where the first sum is over the two different elements of length n. On R 3 — 
{image points} the sequence <Pn converges to a smooth function <P for N — ► 
oo. Convergence follows because at level N the strengths of the new poles are 
suppressed by at least a factor of q N ~ x , where q — supi a »/ (r^ — Oj) < 1. 
Note also that all image poles in 5, lie in fact in the interior of the concentric 
but smaller sphere of radius a- := af/(rij — a,j). Cutting out the interiors of 
S(a' i7 Ci) i = 1,2 thus leaves the spheres Si with small collar neighborhoods 
the two sides of which are isometrically mapped into each other by U (or 7j). 
Using two copies of the manifold so obtained we can pairwise identify these 
collar neighborhoods using these isometries so that an Einstein-Rosen manifold 
with two bridges results. Their topology is that of the twice punctured "handle" 
S 1 x S 2 with each puncture corresponding to an asymptotically flat end. This 
construction generalizes to any number N of holes (or bridges), where as manifold 
one obtains the twice punctured connected sum of N — 1 handles. (For the notion 
of connected sums see e.g. Giulini 1994.) For two holes of equal mass one may 
also just identify Si and S 2 and get Misner's wormhole (Misner 1960) if one 
uses inversions of the first kind, or its non-orientable counterpart if one uses 
inversions of the second kind (Giulini 1990). Both manifolds just have one end. 
In the second case one has the additional possibility to just close "close-off" 
the spheres Si individually by identifying its antipodal points using /j (Giulini 
1992). The manifold has the topology of the once punctured connected sum of 
two real projective spaces RP 3 . It is orientable and has only one asymptotically 
flat end. It can be seen as the generalization to two holes of the once punctured 
single RP 3 obtained above. This construction also generalizes to any number N 
of holes and one obtains the once punctured connected sum of N RP 3 's. These 
manifolds are doubly covered by the iV-bridge manifolds discussed above. 



3.3 Analytic Expressions 

In the case of two holes there exists a geometrically adapted coordinate system 
- so called spherical bi-polar coordinates - which allows to write down explicit 
expressions. We take ai = a 2 = a, ci = de z and c 2 = —de z . Taking the 
Oj's equal means that the holes are of equal size (individual mass). We thus 
consider a two parameter family of configurations labeled e.g. by mass (overall or 
individual) and separation. All image poles are on the ^-axis whose strengths a n 
and locations d n (positively counted z coordinate) satisfy the coupled recursion 
relations 

a i j a 2 
a n = a n -i—— , d n = dT-j—; , (41) 

d + dn-l d+dn-l 
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where the upper (lower) sign is valid for inversions of the first (second) kind. Us- 
ing instead of a, d the parameters c, /Jo defined by a := c/ sinh/_i , d :~ c coth/_<o 
we can solve the recursion relations by 



a n = 



sinh nfio ' 

for the upper sign, and for the lower sign 

c 



d n = c cothn^o, 



a„ = 



sinh n/Uo 
c 

cosh n/xo ' 



d n = c cothn/io 
d n = c tanh n/io 



for n even, 
for n odd. 



(42) 

(43) 
(44) 



In the xz-plane we introduce bi-polar coordinates via exp(/j — 777) = (£+c)/(£ — c) 
with £ = z + ix. By construction the lines of constant /i intersect those of 
constant 77 orthogonally. Both families consist of circles; those in the first family 
are centered on the z-axis with radius c/ sinh/j at \z\ — c coth/j, and those in 
the second family on the x-axis with radius c/ sin 77 at \x\ = c cot 77. Rotating this 
system around the z-axis with azimuthal angle 4> leads to the spherical bi-polar 
coordinates. Explicitly one obtains 



sin 77 cos ( 



V = c- 



sm 77 sin < 



z = c 



sinh /1 



cosh fi — cos 77 cosh \x — cos 77 

Together with (42-44) this gives 

a n [cosh [i — cos 77] x l 2 



cosh fi — cos 77 



||x ± d n e z \\ [cosh(/x ± 2n/i ) - e cost?] 1 / 2 ' 



(45) 



(46) 



where e = 1 if one uses inversions of the first kind and e = — 1 if one uses those 
of the second kind. The final expression for the metric in (/j, 77, 0)-coordinates 
can now be written down: 



71=1 



|x + d„e 2 



llx - d n e z 



-1 4 



^ (cosh(/j + 2n/j ) — e" cos rj) 
.neZ 



-1/2 



dx • dx 



(d/_i 2 + rf?7 2 



sin 2 77 < 



(47) 



2 ). (48) 



It nicely exhibits the isometries (//, 77, 0) 1— > (/j+2/xo, »7, </>) for £ = 1 and (/j, 77, 0) 
(/1 + 2/j , it — r],4>) for e = —1. The extrinsic curvature matrix for the surfaces of 
constant \x with respect to an orthonormal basis in 77 and 4> direction is given by 
2<P~ 3 d<I>/dn times the unit matrix. Hence K has only a trace part (the surfaces 
of constant \x are totally umbillic) and vanishes iff fi = ±fio- Hence in both cases, 
e = ±1, the apparent horizons are also totally geodesic (this we already knew 
for e = 1). 
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Next we turn to the expressions for the masses. We shall follow Lindquist 
(1963) and define the mass of the first hole by appropriately applying (35): We 
sum all the "bare masses" 2a^ in Si, each enhanced by an interaction factor 
1 + Xi which includes the interactions of each pole in Si with any pole in S2, 
but not with poles in Si . This we write as 



(49) 



with the obvious meaning of "g" . Since mi = ni2 we write m for the individual 
mass and M for the overall mass. The latter is just the sum of all 2<Xj. Using 
(42-44) one obtains (quantities referring to e = — 1 carry a tilde) 




m 



00 00 
^—J sinhnun ^— ' smhn/in 

n=l n=l 



for e = 1, and for e = — 1 

In 



m - 



h 2c > ■ (51) 

sinh2nun ^— ' cosh(2n + l)nu ' 

n=l ra=0 x ' 

V 1 X 1 

M = 4cV— — +4cV — . (52) 

^— ' sinh2nun , cosh(2n+ l)n//n 

n=l n=0 x ' 

As mentioned above, we define the distance of the holes as the geodesic 
distance of the apparent horizons /i = ±/Uo- The shortest geodesic connecting 
these two surfaces is n = n. For e = 1 its length, I, may be expressed in closed 
form: 

I = 2c(l + 2mfi ), (53) 

with m from (50). I have not been able to find such a compact expression in the 
case e = —1. 

Like I, many quantities of interest cannot be evaluated in closed form. In 
these cases it may be useful to expand in powers of m/l. Numerical studies show 
that additional outer apparent horizons form (i.e. the holes merge) for values 
above m/l ~ 0.26 (Smarr et al 1976), so that good convergence holds up to the 
merging ratio. 



Comparing e = ltoe=— 1. We have seen that mathematically these two 
cases differ by allowing different topologies. But are there more physical aspects 
in which they differ? A natural question is how for fixed "physical" variables 
m = fh and I — I the total energies M and M differ (Giulini 1990). One finds 
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showing that for e = —1 the holes are slightly tighter bound (i.e. they at- 
tract stronger) , although the additional energy gained until merge is only about 
10 _2 M. This result is qualitatively unchanged if one uses Penrose's instead of 
Lindquist's definition of mass. 

Another difference shows up in the deformation of the apparent horizons upon 
(adiabatic) approach of the two holes. One can define an intrinsic deformation 
parameter as follows: Regard (rj, <f) as polar coordinates. The poles are the zeros 
of the Killing field d^. Define C v as twice their geodesic distance. Among the 
orbits of d,/, is one of greatest length, C<$,. The deformation parameter is D := 
(C n - C^/Cr,. One obtains (Giulini 1990) 



The power of 2 in (56) seems in conflict with the usual "tidal-force" interpre- 
tation. The shapes themselves are also different. Like eggs with the thick ends 
pointing towards each other in the first case, and prolonged symmetrically (with 
respect to reflections on the equator i] = n/2) in the second. 
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